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ABSTRACT 

We provide a comparative study of the Higgs boson mass measurements based on two 
approaches to the dileptonic decay of W bosons produced by the Higgs boson decay, one 
using the kinematic variable M^"*^ and the other using the Mr2-assisted on-shell 
reconstruction of the invisible neutrino momenta. We find that these two approaches can 
determine the Higgs boson mass with a similar accuracy for both of the two main 
production mechanisms of the SM Higgs boson at the LHC, i.e. the gluon-gluon fusion 
and the weak vector boson fusion. We also notice that the Higgs signal distribution for 
the gluon-gluon fusion becomes narrower under the Mt2 cut, while the corresponding 
background distribution becomes flatter, indicating that one might be able to reduce the 
systematic uncertainties of mass measurement with an appropriate Mt2 cut. 
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1 Introduction 



The utmost target of the LHC is to discover the Higgs boson and study its fundamental 
properties to estabhsh the mechanism of electroweak symmetry breaking as the origin of 
particle masses [1]|2] . The unsuccessful search at the LEP experiment set a lower bound on 
the standard model (SM) Higgs mass at 114.4 GeV (95 % CL.) |3], while the analysis of the 
electroweak precision data indicates a relatively light SM Higgs boson with rriH ~ 185 GeV 
at the 95 % confidence level Combined with the recent Tevatron search excluding 158 
GeV < rriH < 175 GeV [5], we anticipate that the SM Higgs boson most probably lies in 
the mass range between 114.4 GeV and 158 GeV. 

At the LHC, the Higgs boson is mainly produced by the gluon-gluon fusion (GGF) 
mechanism, and the second most important source is the weak vector boson fusion (VBF) 
process [6l[7] . While the production cross section of GGF is about 10 times larger than that 
of VBF for the Higgs mass region 114.4 GeV < rriH < 170 GeV, VBF has its advantage 
in a kinematic structure containing two forwarding tagging jets with a large rapidity gap, 
which can be exploited to isolate the Higgs boson signal from backgrounds and to study 
the signal properties [HHTO]. 

The specific search strategy of the Higgs boson at the LHC depends on its mass 
and decay pattern. According to the combined study on the expected Higgs discovery 
significance at ATLAS with 10 fb~^ luminosity [1], when the Higgs boson is lighter than 
130 GeV, the main search channel is GGF with the Higgs boson decaying into two tau 
leptons or into two photons. On the other hand, for 130 GeV < rriH < 150 GeV, the 
following three channels are available with a similar significance: (z) GGF with the Higgs 
boson decaying into two Z bosons which subsequently decay into four charged leptons, 
{a) GGF with the Higgs boson decaying into two W bosons, which subsequently decay 
into two charged leptons and two neutrinos, and {in) VBF with the Higgs boson decaying 
into two W bosons, which subsequently decay into two charged leptons and two neutrinos. 
Finally, for 150 GeV < niH < 190 GeV, Higgs boson decay into two Z bosons is relatively 
suppressed, while channels (ii) and {iii) remain as the dominant search channels [H]. 
Therefore, dileptonic W boson decays play a crucial role in the Higgs boson search at the 
LHC when the Higgs boson weighs between 130 GeV and 190 GeV. Furthermore, if one 
can improve the efficiency of dileptonic channels, they might play an important role even 
for TTiH lighter than 130 GeV. 

In dileptonic W boson decays, there are two invisible neutrinos which make a direct 
reconstruction of the Higgs boson mass impossible. To overcome this difficulty, one can 
consider various kind of transverse mass variables^ A well-known example is the transverse 

*The general concept of transverse mass has been introduced and studied as early as in Ref. [T^]. In our 
work, we concentrate on the recently proposed transverse mass variables such as M!j^^°, Mt2, and rn^^°^ . 
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mass of a boson pair in the process H — )■ WW — )■ Ivl'v' {i, i' = e, fi and u, u' = Ug, v^): 



Ml{WW) = ml, + ml,, +2 [J pf + ml,J p^'^' + ml,, - p« ■ p^/ 
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where mx and p^ denote the invariant mass and the transverse momentum, respectively, 
of X = uu'. Obviously Mt{WW) is bounded by mn, and this upper bound is saturated 
when W and uu' have the same rapidity, rju, = r]^^,, where r] = ^ln{E +pl)/{E — Pl) 
for the longitudinal momentum p^. Therefore, if Mt{WW) is correctly reconstructed 
event by event, the Higgs boson mass can be read off from the endpoint value of Mt{WW) 
distribution. However, there is an obstacle to this approach: m^,y' cannot be experimentally 
determined, although p^"' can be deduced from the missing transverse momentum of the 
event. One way to bypass this difficulty is to simply take m^^/ = mw [13], and consider 
the distribution of 

M— ^MHPyW^)L^^,_^^,, (2) 

which would provide a good approximation to the true value of Mt{WW) if the W boson 
pair were produced at near-threshold. However, in reality, a sizable number of events are 
not close to such a threshold, and as a result, M^^^'^^ is not bounded by m^j anymore. Still, 
the detailed shape of its distribution has a certain correlation with mj^, so M^^''™^ has been 
widely used in the previous studies of the Higgs boson search and mass measurement [UIHl 

m. 



Recently, it has been noticed that an alternative transverse mass variable 

Mr^^Mr(iyiy)L^^,^o, (3) 

may determine the Higgs boson mass more accurately than M^"^"^"^ does, since it is bounded 
by mn and thus could have a stronger correlation with mn [HllIS]. Obviously, for each 
event 

M;^™" < Mt{WW) < mn, (4) 

and the upper bound of M^^'^ is saturated by the events with rju, = rji^y, and mj,,^/ = 0; 
therefore, the endpoint value of M)f"^ distribution indeed corresponds to m^- 

There is another widely discussed transverse mass variable, Mt2, which is defined 
for a generic event with two identical invisible particles [TB] and was applied recently to 
the mass measurement of supersymmetric particles [161 [17]. For the Higgs boson event 
H WW i{p)u{k)i'{qy{l), Mt2 is given by 



Mt2 = niin 

kT+lT=)^>r 



max 



{m«,m(^^}], (5) 

where Adf'^ and M^,^^ are the transverse masses of the two leptonically decaying W bosons: 
[M^'^y = 2 (IpTllkrl - PT ■ kr) , {M^^^f = 2 (|qT||lT| - Qt ■ It) • (6) 
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This Mt2 has an endpoint at mn/'^ when mn < 2mvK [IE], suggesting that Mt2 also might 
be useful for the Higgs boson mass measurement. 

In fact, with Mt2 and additional on-shell constraints, one can approximately recon- 
struct the invisible particle 4- momenta in each event [I9]. This method of "Mr2-assisted- 
on-shell (MAOS) reconstruction" of invisible particle momenta has been applied to the 
Higgs boson event H — WW l{p)v{k)l' {q)v' {I) [18] in order to examine the invariant 
mass variable 

(m^"°^)' = {p + q + fc'"""^ + r'^'^'f , (7) 

where k^^^ and Z™*^"*^ are the reconstructed neutrino 4-momenta. It was then argued that 
the distribution of this MAOS Higgs mass exhibits a peak at the true Higgs boson mass. 
Furthermore, the reconstructed MAOS momenta k™^^ and l"^^^ are closer to the true 
neutrino momenta for the events near the upper endpoint of Mt2- As a result, the rnf[^^ 
distribution can exhibit even a narrow resonance peak at niH when a strong Mt2 cut is 
employed. These observations suggest that imJ^^^ might provide a powerful tool to probe 
the Higgs boson mass through dileptonic W boson decays. 

Since it has been claimed that various kinematic variables are useful for the Higgs 
boson mass measurement, it would be instructive to perform a comparative study of the 
efficiencies of those variables in the LHC environment. In this paper, we wish to perform 
such a comparative study for M^^^, Mt2i and m^^^, and explore the possibility to improve 
the efficiency of the analysis. This paper is organized as follows. In Sec. [2l we discuss 
the distinctive features of M^^'^, Mt2, and rn^^^ for dileptonic W boson decays, mainly 
focusing on the possible correlations between these three kinematic variables. From the 
discussion in Sec. |21 it is evident that Mt2 is less efficient than "'^ as a main observable 
for the Higgs boson search. However, Mt2 might be a useful cut variable in the analysis 
using M^^'^ or m^^^ as the main observable. As for the comparison of M^^^ and m^^^, 
it is not obvious which one is more efficient since each variable has its own virtue and 
weak point. We thus perform in Sees. |3]and|4]a detailed comparative study of M^'^'^ and 
^maos_ Our results show that m'^'^^ and M^^^ will eventually have a similar efficiency for 
both the GGF and VBF Higgs productions at the LHC. In this study, we also examine 
the possible role of Mt2 as a cut variable, and find that the Higgs signal distribution for 
the GGF channel becomes narrower under the Mt2 cut, while the background distribution 
becomes flatter. Although an estimate of any systematic uncertainty is beyond the scope 
of this paper, such behavior of the signal/background distributions under the Mt2 cut 
indicates that one might be able to reduce the systematic errors in the GGF channel with 
an appropriate Mt2 cut. We summarize our conclusions in Sec. |5l 
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2 Distinctive features of the transverse mass variables 



Unlike the Higgs boson decay into two photons or two tau lepton^, the Higgs boson mass 
cannot be reconstructed in the event H — — > ^{p)u{k) due to the missing 

neutrinos. One then often considers the transverse mass of the W pair, Mt{WW), defined 
in Eq. ([1]), as the main observable to probe the Higgs signal event [Il[8l[T3l[Tl]. If one 
could determine the correct value of m^yi event by event, the resulting distribution of 
Mt{WW) would have an upper endpoint at mn- However, 171^,^/ is not available, in general, 
and there have been two proposals to fix the unknown m,^j,/. The first one is to choose 
niui^' = mu' [I1IH1II3], and consider the distribution of 



M— ^MHiyW^)|„,^,^,_^^,. (8) 

This has been motivated by the observation that, for a Higgs boson mass close to 2mw, the 
W bosons are produced at near-threshold and almost at rest in their center-of-mass frame, 
for which m*'j^° = m«/, where m*'""/'^ is the true value of niyyi in the event. Although this 
threshold approximation is not necessarily suitable for generic events and therefore M^"^"^^^ 
is not strictly bounded from above by niH-, still the shape and range of its distribution can 
provide information on the Higgs boson mass. The second proposal [T^ is deduced from 
the simple relation 



M^^^ = Mt{WW)\^^^,=, < MT{WW)U^^,=„^-r < TJlH. (9) 



It was reported that one can obtain a more accurate Higgs boson mass by using M^^^ 
rather than using M!p'^'^°^ p3|[T5]. 

For the event H WW^*'^ i{p)u{k)i'{qy{l), one may exploit the event variable 
Mt2 which has been designed for a general situation with two invisible particles (with the 
same mass) in the final state [16]. The variable is defined as 



Mt2 = niin 

i(.T+iT=t>r 



max {m^^\ M^^]] , (10) 



where M^^^ and M!p^ are the transverse masses of the two leptonically decaying W bosons: 

M^'^y = 2 (IpTllkrI - Pt ■ kr) , {M^^^y = 2 (Iq^lllTl - qT ■ It) • (11) 

For a specific set of events for which the Higgs boson has vanishing transverse momentum, 
]^rp = —pT — qr, one can use the following analytic expression of Mt2- 

= 2(|pT||qr| +PT-qT) (12) 



^ In the di-tau channel, each tau decays to hadron (or lepton) + a neutrino(s). Although the final state 
neutrinos make a full reconstruction of the event impossible, it is known that the di-tau invariant mass can 
be determined well in collinear approximation |20| . 
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to find the upper bound 



< (|PT| + |qT|)' < (13) 



Although the above relation is derived for the specific event set with vanishing Higgs 
transverse momentum, it is likely that the same upper bound applies for generic Higgs 
boson events with nonvanishing Higgs boson transverse momentum. We found through 
numerical analysis that it is indeed the case and the inequality Mt2 < / 2 remains true 
for generic dileptonic Higgs boson events H — )■ H^l^^*'' — )■ iv^v' . [See the later discussion 
of the inequality between Mt2 and M-JT"" in Eq. f l2T|) .] On the other hand, in the case that 
both W bosons are on shell, Mt2 is bounded by m-w. Combining this with (|T3|) . one may 
find that the M^^ of the Higgs boson decay is bounded by 



j^^max ^ 1 for mu < 2mw 

\ rriw for niH > 2mw, 

implying that Mt2 can provide information on the Higgs boson mass through its endpoint. 

Apart from M^''^™^, M^^^, and Mj'2 in dileptonic W boson decays, there has been 
a proposal to reconstruct an invariant mass by using the MAOS reconstruction of the two 
neutrino momenta [18]. In MAOS reconstruction, the transverse momenta are defined as 
the ones that determine the value of Mt2, i-e. 

Mt2{Pt, qr, I^t) = M^'\pt, = M^'\qT, 1^°' = j^t - k^"'), (15) 

while the longitudinal and energy components are obtained from the following constraint^: 

It was noticed that the MAOS 4-momenta determined as the solutions of fllSp and fll6l) 
provide a reasonable approximation to the true neutrino 4-momenta, and they become 
closer to the true momenta for the near-endpoint events in the Mt2 distribution [T8],[T9]. 
Then the distribution of the following invariant mass (= MAOS Higgs boson mass), 

(m^"°^)' = {p + q + A;'"""^ + r^^'f (17) 

exhibits a peak at the true Higgs boson mass, and the peak shape becomes narrower if 
one imposes an event cut selecting the events near the endpoint of Mt2- Although such 
an Mt2 cut can cause the rnJ^^"^ distribution to have a resonance peak at ttlh and also 
might eliminate some of the backgrounds, it can also sacrifice the number of the Higgs 



^Note that one could use mw instead of Mt2 in the second constraints if both of the W bosons are on 
shch. Our scheme to use Mt2 can be apphcable irrespective of whether the W bosons are on sheU or off 
shcU. 
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signal, thus worsening the statistical significance of the mass measurement. Therefore for 
a given luminosity a careful analysis is required to see if the efficiency of the Higgs mass 
measurement can be improved by a proper cut. 

To proceed, let us consider a specific set of events for which the W boson pair has 
a small transverse momentum: pj!^^ ^ 0. Although it does not cover the full event set, 
this subset of events reveal some of the essential features of the kinematic variables under 
discussion. For those events with small , we have 

P^T = -(pt + qr), 

and then the following simple expressions of Mt2 and MAOS momenta are available [T8UT9] : 



kmaos 
T 



imaos 



qT 



jj^maos 



„ rmaos 

-Pt , II 



kmaos I 
I T I 

IptI 

I imaos I 

Iqrl 



PL 



= 2 (IpTHqrl + Pt ■ qr) 



With these expressions of Mt2 and MAOS momenta, it is straightforward to find 



(m 



H 



Pt 



V qT 

Ipt + qTl 



qT 



Pt 



2 (IPTHqrl + Pt ■ qT) + m^w) , 



Pt + qT|^ + m%. 



(19) 



Since 



nil, 



2 (|pT||qT| + Pt ■ qT) < |pT + qT| 

2(|p||q|-p-q) > 0, 



(20) 



both m^^°^ and M^j^^'^ are bounded below by 2Mt2, and therefore 

2Mt2 < M^""™ < tuh , 2Mt2 < m^^°^ . 



(21) 



Although the above inequalities between Mt2, M^™ and m^^^ could be derived 
analytically only for the events with p^'^ = 0, we find numerically that they hold true 
also for the events with p^'^ 7^ 0. In Fig. [1] (left panel), we show the distributions of these 
three kinematic variables for the true Higgs mass mn = 160 GeV. Here we use the Higgs 
boson events produced by the GGF process at the LHC, while incorporating the detector 
effects with the fast detector simulation program PGS4 [21]. Figure [T] (left panel) shows, 
first of all, that both M^'^^ and 2Mt2 have a common upper endpoint at ruH = 160 GeV, as 
we have anticipated in the above discussion. However, the M^'^^ distribution is significantly 
narrower than the 2Mt2 distribution, as implied by the first inequality in Eq. ( 12T|) . This 
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indicates that Mjf"'^ is more sensitive to mn than Mt2, and thus is a more efficient variable 
to determine rriH- If one compares M^'^'^ to m^^°^, Fig. [T](left panel) shows that m^^^ has 
a peak at m^^, while M^''™ has an endpoint at mu- With this feature, one might expect 
that m^'^^ is more sensitive to mH than M^^° . However the M^™^ distribution around niH 
is significantly narrower than the m^^^ distribution, and thus one needs a detailed analysis 
to see which of M^™ and rn^^^ is more efficient for the determination of mu- In the next 
two sections, we perform such an analysis for both the GGF Higgs production (Sec. [3]) and 
the VBF Higgs production (Sec. HD at the LHC. 



3 Measuring the Higgs boson mass in GGF 

In this section, we consider the SM Higgs boson production at the LHC through the GGF 
process and its decay into two leptonically decaying W bosons. To investigate the experi- 
mental performance of m^*^"*^ and M^''™, we use the PYTHIA6.4 event generator [22] at the 
proton-proton center-of-mass energy of 14 TeV. The generated events have been further 
processed through the fast detector simulation program PGS4 [21] to incorporate the de- 
tector effects with reasonable efficiencies and fake rates. The dominant background comes 
from the continuum qq , gg ^ WW — )■ /z//V process, and we also include the tt background 
in which the two top quarks decay into a pair of W bosons and two b jets. For the details 
of the MC event samples of the SM Higgs boson signal and the two main backgrounds, and 
also the employed event selection cuts, we refer to Ref. [TH] . 

In Fig. in we show the m^^°^ and M^^^ distributions of the signal events for mn = 160 
GeV without (left panel) and with (right panel) the Mt2 cut: Mt2 > {,Mt2T''^ = 60 GeV. 
We use the same event samples at the detector level as in Ref. [18] but do not apply the 
event selection cuts other than the Mt2 cut. First we note that the M^'^'^ distribution is 
bounded above by rriH = 160 GeV, while the m^^°^ distribution is peaked around it, as 
noticed in the discussion of the previous section. Under an Mt2 cut, both distributions are 
bounded from below by 2 (Mt2Y^^ = 120 GeV with a rather good approximation, which 
can be understood by the inequalities in (^T^. Incidentally, in the m^^^ distribution, the 
region above mu is less sensitive to the Mt2 cut than the region below mn- We also see 
clearly that the M^^'^ distribution is narrower than the m^^^ distribution as implied by 
Eq. 

Next, we proceed to perform the template likelihood fitting to the M^^'^ and m^^°^ 
distributions. Here, a template means a simulated distribution with a trial Higgs mass 
which, in general, is different from the nominal one used to generate the data. For each 
pseudoexperiment distribution with the nominal Higgs mass mn, 11 templates are gener- 
ated with the trial Higgs masses between rriH — 10 GeV and rriH + 10 GeV, in steps of 2 
GeV. For example, in each frame of Fig. |21 the solid line shows the template distribution 
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Figure 1: The signal m^^^ (black solid line) and M^^^ (red dashed line) distributions at 
the detector- level without (left panel) and with (right panel) the Mt2 cut, Mt2 > 60 GeV. 
The Higgs boson mass is taken to be 160 GeV, and we use the same event samples as in 
Ref. [ig. 

when the trial Higgs mass is the same as the nominal one. Each template is normalized 
to the corresponding pseudoexperiment distribution. For the definitions of the likelihood 
between a pseudoexperiment data distribution and a template, and also for more on the 
Higgs masses and l-cx errors obtained by fitting the log likelihood distributions, we again 
refer to Ref. [18]. 

Empirically, the template likelihood fitting, which is used to estimate the efficiency of 
the Higgs mass determination in this work, provides a better result when the width of the 
distribution is narrower. On the other hand, compared to the peak of the distribution which 
is generally easier to determine, the endpoints are more vulnerable to detector smearing, 
backgrounds and low statistics. Even though the peak of the M^^^ distribution does not 
have a strong correlation with the input Higgs mass like that of the m^^^ distribution, 
the overall shape has a definite correlation with the Higgs boson mass and, moreover, its 
narrower width might result in a comparable sensitivity to mlff"^. The problem is the shape 
of the background distribution which could smear the edge structure of M}^^'^. In Fig. HI 
we have shown the m^^°^ (left panels) and M^^^ (right panels) distributions for the Higgs 
signal and the backgrounds. We use the same detector-level event samples as in Ref. [18], 
and all the event selection cuts are applied. 

Our MC study shows that, generically, the background distributions for both M^'^^ 
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Figure 2: The m^'^^ (left frames) and the M^^'^ (right frames) distributions at 10 fb~^ 
luminosity. In each frame, the shaded region represents the backgrounds {ti and W^VT) 
and the event selection cuts without (upper frames) and with (lower frames) the Mt2 cut 
{Mt2 > 66 GeV) imposed. The Higgs boson mass is taken to be 160 GeV, and we again 
refer to Ref. [18] for the details. 



and rriff^^ become flatter if we impose a stronger Mt2 cut. On the other hand, a stronger 
Mt2 cut causes the signal distributions to have a narrower shape (see Fig. [1]). This behavior 
of the signal and background under the Mt2 cut suggests that one might be able to improve 
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Figure 3: The bands showing the 1-a deviation of the statistical error for the Higgs boson 
mass determined by the m^^"^ and M^™*^ distributions in the GGF process gg ^ H ^ 
lul'hi'. The dots and hues denote the Higgs boson mass obtained by the hke- 
hhood fit to the distributions. We use the same event sample and the cut procedure as 
in Ref. [18], where the right (left) panel is the result with (without) an optimal Mt2 cut, 
chosen differently for different values of mn- 

the efficiency of the Higgs mass measurement with an appropriate Mt2 cut, particularly 
reduce the systematic uncertainties associated with various origins such as the detector 
energy resolution, fit routines and poorly estimated backgrounds. As for statistical error, 
one needs a detailed analysis, as the improvement due to better shapes of signal and 
background distributions can be compensated by the reduced statistics. Indeed, for an 
integrated luminosity ~ 10 fb~^, we find (Fig. E]) that there is no appreciable improvement 
of the statistical error gained by the Mt2 cut, although the situation might be different at 
higher luminosity. 

In Fig. 131 we show the Higgs boson masses and the l-a (statistical) errors obtained 
by the likelihood fit to the m^^^ and M^^'^ distributions in the GGF process gg ^ H ^ 
WW^*^ Ivl'v'. We see that the reconstructed mass is rather close to the input one. For 
the sake of comparison, we have not employed an Mt2 cut in the left frame since the cut was 
not included in the original suggestion of the M^^° variable [H] ; however, it is employed 
in the right frame. We observe that the efficiency of m^^°^ is comparable to that of My "'^ 
up to systematic errors, which are not considered in our analysis. In particular, when the 
Mt2 cut is employed, it is hard to say which one shows a better efficiency. Still, we note 
that m^^°^ could be a better choice when the distribution is spoiled by some unknown 
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backgrounds beyond the SM, since the peak is less vulnerable to unknown backgrounds 
than the endpoint. 

4 Measuring the Higgs boson mass in VBF 

In this section, we turn to the VBF process, which is the second most important production 
channel of the SM Higgs boson at the LHC. The characteristic feature of the process is 
the existence of the two forward tagging jets with suppressed hadronic activity between 
them and the central Higgs decay products. Thanks to the exchanges of the colorless 
vector bosons in the process, these features lead to a fairly clean environment with well- 
isolated signal events in a low background. Therefore, the VBF process is useful to measure 
properties of the Higgs boson in a hadron collider environment. 

The same variables introduced in the GGF process could be used to determine the 
Higgs boson mass in VBF, since the tagging jets provide additional information without 
touching the decay of the Higgs boson itself. To investigate the experimental performance 
of the MAOS Higgs mass at the LHC through the VBF process, we have generated MC 
event samples of the signal H + 2j events and the main ti background events by PYTHIA6 . 4 
[22] . The subleading background, i.e. electroweak WW + 2j event showing a similar 
characteristic, has been generated using MadGraph/MadEvent [23] • The generated events 
have been further processed with the fast detector simulation program PGS4 [21], which 
approximates an ATLAS or CMS-like detector. The PGS4 program uses a cone algorithm 
for jet reconstruction, with a default value of the cone size Ai? = 0.5, where Ai? is a 
separation in the azimuthal angle and pseudorapidity plan^. We note that the 6-jet tagging 
efficiency is introduced as a function of the jet transverse energy and pseudorapidity, with 
a typical value of about 50 % in the central region for the high energy jets. 

Proceeding in a similar way as in Ref. [8], we have imposed the following basic selection 
cuts on the Higgs signal and the backgrounds: 

• Preselection cuts 

• 6-jet veto 

• Tag-jet conditions (forward jet tagging, leptons between jets, central jet veto) 

• 7*/Z+ jets, Z — )■ rr rejection 

^The pseudorapidity is defined as i] = — lntan(0/2) for the angle 9 between the particle momentum 
and the beam axis. Note that the pseudorapidity for a massless particle is equal to the rapidity defined as 

in = \HE + pL)l{E~PL). 
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Since the final state consists of two fiard jets, two cliarged leptons and significant missing 
transverse energy, the preselection cuts should be the minimal conditions to be imposed. 
Specifically, we require at least two jets with > 20 GeV and |?7| < 4.8, and two isolated, 
opposite-sign leptons (e or /i) with px > 15 GeV, and Ii^tI > 30 GeV. The fe-jet veto is 
used to exclude the main tt background which contains two b jets. The tag-jet conditions 
should reflect the characteristic features of the VBF process: two forward tagging jets 
with suppressed hadronic activity between them and central Higgs decay products. In this 
work, we define the tagging jets as the two highest pt jets in the event while rejecting the 
event if they are in the same hemisphere. Explicitly, we have imposed the following tag-jet 
conditions. 

— Forward jet tagging: the two hardest jets are required to be in opposite hemispheres 
(?7j^ X ?7j2 < 0) and to have a pseudorapidity separation lArjj-^j^l > 3. 

— Leptons between jets: both leptons are required to be between the tag jets in pseu- 
dorapidity. 

— Central jet veto: the event is rejected if it contains any extra jet (in addition to the 
two forward tagging jets) with p^ > 20 GeV in the pseudorapidity region \ri\ < 3.2. 

The last basic selection cut is to reject the Drell-Yan background involving the variable 
[see Eq. (1)]: 

m'^" = MT{WW)U=m.^=o. (22) 

Though the Drell-Yan background has been neglected in this work, we have included the 
rejection cut to see its effect on the signal and the other background events under con- 
sideration. For the Z — )■ rr rejection, the di-tau invariant mass is constructed using the 
collinear approximation, rejecting the events with \MrT — fnz\ < 25 GeV. 

It is well known that the spin-zero nature of the Higgs boson makes two charged 
leptons come out in the same direction, making its opening angle in the Higgs signal 
smaller than that in the background; see Fig. HI To incorporate this spin information, in 
addition to the basic selections, we have further required 

A$„ < A<Dr , Arju < Ar^T , (23) 

where A$;/ and Arju are the transverse opening angle and pseudorapidity gap between the 
charged leptons. 

In the literature, the dilepton invariant mass Mu and the dijet invariant mass Mjj for 
two tagging jets are also used for the event cut. In Fig. [3, we show the Mu and Mjj distri- 
butions for both the Higgs signal and the tt background. As in the GGF process discussed 
in the previous section, one might employ Mj'2 as a cut variable in the VBF process also. 
In the Higgs signal with the two tagging jets and also in the main ti background with two 
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Figure 4: Transverse opening angle A^u (left panel) and pseudorapidity gap Arjn (right 
panel) between the charged leptons in the VBF process (m^ = 160 GeV) and the back- 
ground events. The preselection cuts are used for the event selection. 




Figure 5: The Mu (left panel) and Mjj (right panel) distributions of the signal and the tt 
background, after applying the preselection cuts. niH = 160 GeV is taken for the signal. 
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Figure 6: The (left panel) and (right panel) distributions of the VBF process 

{niH = 160 GeV) and the ti background, after applying the preselection cuts. 



b jets, one can construct two sorts of — one using the transverse masses of the ly-pair 
system and the other using the VT-pair + two jets system: 







= min 


max 








= min 

]i.T+iT=t>r 


max 




where M^^^^^ 


and M^""('^ 


are the transverse 


masses 



sub (2) ~ 



Ml 



(24) 



respectively. The variable is defined without any combinatorial ambiguity 
and bounded from above by min(mvK, fnu/'^) [see Eq. (|T^ ] in the signal and by mw in 
the background. On the other hand, there is a combinatorial uncertainty in constructing 
M^2 since, in each selected event, there are two possibilities in associating two jets with 
two charged leptons. Between the two possibilities, we choose the combination giving the 
smaller M^2- Then, for the ti background, is bounded from above by mt at parton 
level, which suggests that M^2 can be used as a cut variable to eliminate the ti backgrounds 
by requiring M^2 > ^t- For the Higgs signal events, tagging jets tend to have a large M, 



which has been used to select the signals by requiring Mjj > MjJ^. In Fig. [6l we show the 
(left panel) and M^2 (right panel) distributions for the signal and the ti background 
at detector level. In the right panel of Fig. El we observe that the M^2 distribution for 
ti backgrounds has a non-negligible tail above mt, which arises from the fact that one 
(or both) of the tagging dijets used for the construction of M^2 is not the b jet from the 
top quark decay, but a hard jet coming from the initial state radiation or an erroneously 
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Table 1: Cut flows (in fb) for tuh = 160 GeV. 



Selection cuts 


VBF H WW 


tt 


EW WW+ jets 


S/B 


Preselection 
6-jet veto 

Forward jet tagging 
Leptons between jets 
Central jet veto 

m^jy > 30 GeV, Z ^ tt rej. 

Mil < 85 GeV 

Mjj > 550 GeV 


33.8 
32.8 
22.5 
20.9 
16.8 
15.1 
14.9 
9.3 


10910.0 
5193.1 
542.1 
347.3 
166.9 
138.5 
62.0 
8.0 


15.9 
15.5 
6.3 
5.6 
3.5 
2.7 
0.8 
0.7 


0.0031 

0.0063 

0.041 

0.059 

0.099 

0.11 

0.24 

1.07 


A$z; < 1.6, Ar]ii < 1.5 


8.4 


5.5 


0.5 


1.39 


M^f > 160 GeV 
M^f > 60 GeV 
A^u < 1.6, Ar]ii < 1.5 


8.9 
3.6 
3.6 


7.1 
1.7 
1.6 


0.7 
0.3 
0.2 


1.14 
1.80 
2.00 


Table 2: Cut flows (in fb) for mn = 140 GeV. 


Selection cuts 


VBF H -> WW 


tt 


EW WW+ jets 


S/B 


Preselection 
6-jet veto 

Forward jet tagging 
Leptons between jets 
Central jet veto 

m^^ > 30 GeV, Z ^ tt rej. 

Mil < 85 GeV 

Mjj > 550 GeV 


17.6 
17.3 
12.2 
11.6 
9.8 
8.5 
8.5 
5.7 


10910.0 
5193.1 
542.1 
347.3 
166.9 
138.5 
62.0 
8.0 


15.9 
15.5 
6.3 
5.6 
3.5 
2.7 
0.8 
0.7 


0.0016 

0.0033 

0.022 

0.033 

0.058 

0.060 

0.14 

0.66 


A$,z < 1.7, Ar]ii < 1.6 


5.4 


5.5 


0.5 


0.89 


M^^^ > 140 GeV 
M^f > 45 GeV 
M^f < 72 GeV 
A$„ < 1.7, Ar]ii < 1.6 


5.5 
2.1 
2.1 
2.1 


7.5 
3.7 
3.1 
2.6 


0.7 
0.5 
0.3 
0.3 


0.67 
0.50 
0.62 
0.72 



reconstructed jet. With this, comparing to Mjj in Fig. [5l we see that M^2 can be as 
efficient as Mjj in discriminating the Higgs signal from the tt background. On the other 
hand. Figs. |6] and [5] show that both the Higgs signal and the ti background have a similar 
shape of M^ distribution, so M^ is not so useful as a cut variable. 
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Table 3: The Higgs boson masses and the 1-a statistical errors obtained by the likelihood 
fit to the 771^^°^ distributions in the VBF processes at 10 fb~^. 



tuh (GeV) 


130 


140 


150 


160 


170 180 190 


Fitted value (GeV) 


132.5 


140.7 


150.8 


161.0 


170.6 183.4 188.6 


1-0" error (GeV) 


7.5 


7.0 


4.8 


3.7 


3.8 6.2 6.8 



With the above observations, in our event selection we have applied Mu and Mjj cuts 
after the basic selection cuts: 

Mu < M[^\ Mjj > M™*. (25) 

Then, more refined Mt2 cuts have been introduced to reduce the number of background 
events further. In Tables [1] and [2l we show how the cross sections of the signal and 
backgrounds change under each selection cut for rriH = 160 GeV and 140 GeV, respectively. 
We note that one may also use the upper cut on the when rriH < 2myy] see the third set 
of cuts in Table O This is because the signal distribution is bounded above by mi^/2 < rriw 
in this case. This upper cut cannot be applied when mn > 2miy, as the distributions 
of both the signal and the background have an equal maximum, mw. 

In Fig. [71 we show the imJ^^^ distributions for the two nominal values of the Higgs 
boson mass: = 140 and 160 GeV. Here (m?}^°^)^ = (pi + A;f^^°'^ + + for 

1^(1)1^(2) —i. /(i)(p^)z/(^)(/c^)/(2)(^p2)z/(^)(/i;2), where the transverse components of /c™^°** {% = 
1,2) are those used to determine M^2-i ^i-nd the longitudinal and energy components are 
fixed by the constraints (pi + fcj"^"'^)^ = (p2 + k^f'^^f = {M^ff and {k'^''°'f = {k^^^^f = 
0. We observe that the rn^^^ distribution has a clear peak at the true (nominal) Higgs 
massQ With the m^^°^ distribution constructed as above, we performed a template fitting 
to determine the Higgs boson mass as in Sec. |3l the GGF case. The estimated Higgs 
masses, together with the l-cx errors obtained by fitting the log likelihood distributions for 
various Higgs masses, are listed in Table El The 1-a deviated value is defined as the one 
that increases — ln£ by 1/2 [21]. We again see that the reconstructed mass is quite close 
to the input Higgs boson mass. 

For a comparison of the efficiency of m^^^ to that of M^^^, we also perform a likeli- 
hood fit of the M^™*^ distribution with the same event sets, including the backgrounds. In 
Fig. in we show the M^™'^ template distributions (solid lines) when the trial Higgs mass is 
140 GeV (left panel) and 160 GeV (right panel). The dots with error bars represent the 
M^^^ distributions for the pseudoexperiment data with the same nominal Higgs masses. 

^The background bump appearing around rri^^"^ > 200 GeV is an accidental consequence of our analysis 
for the VBF case. We observe that it becomes less significant like as in the GGF case for a different choice 
of the bin size. 
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Figure 7: The m^^°^ distributions for the pseudoexperiment data (dots) and the template 
(sohd hue) for = 140 (left) and 160 GeV (right) at 10 fb'^The hghtly and thickly 
shaded regions represent the ti and WW + 2j backgrounds, respectively. 




Figure 8: The M/)f"° distributions for the pseudoexperiment data (dots) and the template 
(solid line) for mn = 140 (left) and 160 GeV (right) at 10 fb~^ luminosity. The lightly and 
thickly shaded regions represent the ti and WW + 2j backgrounds, respectively. 
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Figure 9: (Left panel) The band showing the 1-a deviation error for the Higgs boson 
mass determined by the m^'^^ and M^^'^ distributions in the VBF process qq — )■ qqH with 
H — )■ — > lul'u'. The dots and hues denote the Higgs boson mass obtained by the 
hkehhood fit to the pseudoexperiment distributions. (Right panel) The same as in the left 
frame but without the Mt2 cuts. 

We show our results on the Higgs boson masses and the 1-a errors in the left frame of 
Fig. ini these results are obtained using the likelihood fit to the m]^^°^ and M^''"° distribu- 
tions in the VBF process qq — qqH with H — j- VTl^^*) — lul'u'. We observe that the 
efficiencies of the two methods are comparable to each other, and both of them provide a 
good accuracy for the determination of the Higgs mass. Furthermore, to see the role of 
the Mj'2 cuts in the measurement, we repeat the likelihood fitting of the distributions after 
applying the same event selection cuts except the Mt2 cuts; see the right frame of Fig. [HI 
We find that the final result does not change much since the purity of the event samples 
is already high enough before imposing the Mt2 cuts and also the number of signals is 
relatively smaller than that of the GGF process. Although M^^^ and m^^^ have a similar 
efficiency for the VBF process, we again note that in a situation where there are unknown 
backgrounds beyond the SM which might smear the endpoint of the M}^^'^ distribution, the 
^maos (3Q^i^ ]-,g better choice since the peak is less vulnerable to unknown backgrounds 
than the endpoint. 
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5 Conclusions 



We performed a comparative study of the Higgs boson mass measurements based on two 
kinematic observables, M^^'^, and dileptonic decays of a boson pair. For this, 

we first discussed some features of M^^^ and Mt2, and also of the MAOS reconstruction 
of the neutrino momenta in dileptonic W boson decays and the associated invariant mass 
observable m^^^. It is found that M^'^'^ and m^^^ can determine the Higgs boson mass 
in the range 130 GeV < mn < 190 GeV with a similar accuracy for both of the two main 
production mechanisms of the SM Higgs boson at the LHC, i.e. GGF and VBF, up to 
systematic errors which are not considered in our analysis. Still, it should be noted that 
the m^'^^ distribution has a peak at rriH, while the M^^^ distribution has an endpoint, and 
thus m^^°^ can be a better observable when there are some additional backgrounds due to 
new physics beyond the SM since it is likely that the peak is less vulnerable to unknown 
backgrounds than the endpoint. One might consider an approach using both m^^°^ and 
M^^*^ together to extract maximal information on the Higgs boson mass from experimental 
data. However, our results in Sec. [2] imply that these two variables have a rather strong 
correlation, and therefore it is not likely that such an approach significantly improves the 
accuracy of the Higgs mass measurement. 

Our study suggests that Mt2 can be a useful cut variable in the GGF process, as the 
signal distribution becomes narrower while the background distribution becomes flatter 
under an Mt2 cut. Although systematic errors are not considered in our analysis, such be- 
havior of the signal and background distributions might be useful, particularly for reducing 
various systematic uncertainties in the real analysis of experimental data. In our analysis 
of the VBF process, Mt2 is not a particularly useful cut variable compared to others such 
as the dijet and dilepton invariant masses. 
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